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I. INTRODUCTION 

The theoretical results of atomic parity non- 
conservation (PNC) when combined with the experimen- 
tal results is an important probe of physics beyond the 
standard model of particle physics [1]. There are two 
sources of PNC in atoms, nuclear spin-independent (NSI) 
and nuclear spin-dependent (NSD). The NSI-PNC is well 
studied and experimentally observed in several atoms. 
The most precise measurement till date is in the case of 
atomic Cs 2] ■ The same experiment also indicated a sig- 
nature of NSD-PNC effects. The most dominant source 
of which is the nuclear anapole moment (NAM) , a parity 
odd nuclear electromagnetic moment arising from parity 
violating interaction within the nucleus [3H5]- However, 
there are two other contributions to NSD-PNC, these are 
the NSD electron-nucleus Z exchange interaction and the 
combined effect of hyperfme interaction and NSI electron- 
nucleus Z exchange interaction. 

The parameters describing nucleon-nucleon coupling, 
effect of NAM is subsumed into it, extracted from the 
Cs PNC experiment do not concur with the nuclear data 
6 . This certainly calls for the further investigation of 
the NSD-PNC effects in other atomic systems as well. 
An example of an alternative experiment is the proposal 
to measure the PNC in Ba + ion, suggested by Fortson 
[7J and is in progress at Seattle [5] H] . This experiment 
could lead to an unambiguous observation of NAM in the 
6s 2 Si/2 — 5d 2 D 5 / 2 transition, as the NSI-PNC alone 
does not contribute to this transition. It is important to 
note that the major difficulty to a clear observation of 
NAM is the large NSI signal, which overwhelms the NSD 
signature. The Ra + ion has also been suggested and is 
considered to be an important candidate for the PNC 
measurement [TUl [UJ- Apart from Ba + and Ra + ions 
which are one- valence systems the other promising candi- 
date for PNC, the NAM in particular, measurement is the 
atomic Yb. An enhanced effect of PNC has already been 
reported QJ [T3] in neutral Yb, the 6s 2 1 S - 6s5d 3 D 2 
transition, and for further refinement of the experiment 
is in progress at Berkeley. The 6s 2 <Si/ 2 — 5d 2 D 3 / 2 
transition in Yb 4 ", has also been suggested to reveal the 
NAM signature and is being investigated at Los Alamos 

mug. 



The atomic theory results using reliable and accurate 
many-body methods are key to estimate the expected 
value of PNC transition amplitudes and extracting NAM. 
For the theoretical calculations, the relativistic coupled- 
cluster (RCC) theory [TBI [17] can be of great significance, 
as it is one of the most reliable many-body theory to 
incorporate electron correlation in atomic calculations. 
The RCC has been used extensively in atomic structure 
calculations [10, 18 23J of properties like transition ener- 
gies, hyperfme structure constants, electromagnetic tran- 
sition amplitudes, intrinsic electric dipole moment and 
PNC in atoms. Apart from atomic physics, it has also 
been used with great success in nuclear |24j . molecular 
[2"5] and the condensed matter [26! physics. 

In this work, we employ perturbed relativistic coupled- 
cluster (PRCC) theory to calculate NSI and NSD-PNC 
amplitudes of the [4/ 14 ]6s 2 S 1/2 - [4/ 14 ]5d 2 D 3/2 tran- 
sition in the case of 171 Yb + ion. This is timely as there 
are few theoretical results, Sahoo et al [37] and Dzuba et 
al [28] for NSI-PNC and Dzuba et al [28] and Porsev et 
al [H] for NSD-PNC are the previous works. The NSI- 
PNC results from Ref. [57] calculated using RCC method 
differ substantially from Ref. [25] where the correlation- 
potential-method with sum-over-state approach is em- 
ployed to calculate NSI and NSD-PNC. The NSD-PNC 
results reported in Ref. [5S] are based on RPA and, in 
general, is in agreement with the results reported in Ref. 
|28j . However, the later is based on the sum-over-state 
approach, at the level of PNC matrix elements. The 
PRCC method [3"0r[3"2] employed in present work is dif- 
ferent from the sum-over-states approach. It accounts 
for the all singly and doubly excited intermediate states. 
There are two sets of the cluster amplitudes in the PRCC, 
and the summation over states in the first order time- 
independent perturbation is incorporated in one set of 
the cluster amplitudes. 

The paper is organized as follows. In Section. [TTJ we 
provide a brief description of the theoretical methods. 
The unperturbed RCC equations for close-shell and one- 
valence systems are given to serve as a easy reference. 
The perturbed RCC is then discussed in detail and PRCC 
equations are derived. The expression for E1PNC using 
PRCC wave function and some leading order diagrams 
are also discussed. Results from the work and uncertainty 



estimates are presented and discussed in Section. Ill 



II. THEORETICAL METHODS 
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In absence of PNC interaction the atomic states are of 
definite parity, and we consider these as the eigen states 
of the no-virtual-pair Dirac-Coulomb Hamiltonian [33 
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where oli and /3 are the Dirac matrices, p is the linear 
momentum, Vjv(r) is the nuclear Coulomb potential and 
the last term is the electron-electron Coulomb interac- 
tions. The operator A + projects on the positive energy 
eigenstates to avoid the negative energy continuum so- 
lutions. The Hamiltonian H BC satisfies the eigen value 
equation 

H I}C \^ V )=E V \^ V ), (2) 

where is the exact atomic state of the one- valence 
system and E v is the corresponding energy. Here after, 
for compact notation, we use H to represent H DC . In 
the present work, we use RCC theory with the single 
and doubles (CCSD) excitation approximation to solve 
Eq. ([2]). In RCC, \^ v ) is expressed in terms of the 
closed-shell and one-valence cluster operators, and 
respectively, as 



I*, 
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1 + 



1$, 
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where superscript (0) represents the unperturbed RCC 
operators. The one-valence Dirac-Fock (DF) reference 
state |$„) is obtained by adding an electron to the closed- 
shell reference state, = aJ|<&o). In the CCSD ap- 
proximation, T(°) = if* + T 2 (0) and S<® = + sf. 
Using the second quantized representation 

Ti = ^2tyla a , and T 2 = i ^ ^fa^aja 6 a , (4a) 

a,P a,b,p,q 

S\ = ^ s p v ala v , and #2 = E s p v q a (J p a\a a a v . (4b) 

Here, t.\\ and s\.\ are the cluster amplitudes. The in- 
dexes abc. . . (pqr . . .) represent core (virtual) states and 
vwx . . . represent valence states. The operators T\ [S\ 
) and T 2 (1S2) give single and double replacements af- 
ter operating on the closed(open)-shell reference states. 
The diagrammatic representation of these operators are 
shown in Fig. [l] 

The open-shell cluster operators are then the solutions 
of nonlinear equations [34 

(^ v \Hn + {hVs {0} }\<P v ) = E^(^\s[ 0) \9 v ), (5a) 



FIG. 1. Diagrammatic representation of the single and double 
excitation unperturbed cluster operators in closed shell and 
one- valence sectors. 



where Hys = e _Tt0) 7JNe T<0) is the similarity transformed 
Hamiltonian, = H — ($q\H\Qq) is the normal ordered 
Hamiltonian and E^ u is the attachment energy of the 
valence electron. The operators are the solutions of 
a similar set of nonlinear coupled equations 



(^|ff N |$ >=0, 

«| J ffN|<f>0>=0. 



(6a) 
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The details on the derivation of these equations are given 
in our previous work |35j . 

In presence of PNC interaction atomic states mix with 
the opposite parity states and the total atomic Hamilto- 



ff a = H 



DC 



PNC- 
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Here, A is the perturbation parameter and Hpnc rep- 
resents the any general PNC interaction Hamiltonian. 
It has two components, the NSI and NSD interaction. 
These are 
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where, Gf(= 2.22 x 10~ 14 a.u.) is the Fermi coupling 
constant, Qw and (J>w are respectively the weak nuclear 
charge and the weak nuclear moment of the nucleus ex- 
pressed in terms of neutron and proton numbers, a and 
75 are the Dirac matrices, pn(^) is the normalized nu- 
clear density and / is the nuclear spin. Compared to the 
NSI-PNC, the NSD-PNC require two important consid- 
erations because of the nuclear spin operator I. First, the 
cluster operators in the electron space are rank one opera- 
tors, and second, the atomic states in the one- valence sec- 
tor are eigenstates of total angular momentum F = 1 + J. 

Similar to the unperturbed eigen value equation, Eq. 
([2]), we may write the perturbed eigenvalue equation, sat- 
isfied by the total atomic Hamiltonian, as 



(9) 



where \^ v ) is the perturbed atomic state and E v is the 
corresponding energy. To the first-order in A, \ty v ) = 
+ A|**) and E v = E v + \E*, where the bar in |f *) 
denotes it's parity is opposite to \^ v )- From here on, to 
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FIG. 2. Diagrammatic representation of the single and dou- 
ble excitation NSD-perturbed cluster operators in closed-shell 
and one-valence sectors. The extra line in the T% and 
is to indicate the multipole structure of the operators. 

derive the PRCC equations we consider the NSD-PNC 
into 
Eq. 



interaction Hamiltonian. Using Eq. ( 8b ) , we can rewrite 
© as 



H 
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Here, H: 



NSD 



(G-pn' w /V2) J2i a ipN( r ) is the electronic 
part of Hp^Q. While writing above equation we have 



I -"PNC I *i 



0, as -ffpNc i s an °dd parity 



used El = 

operator it connects opposite parity states only. In the 
PRCC theory, the perturbed wave function is expressed 
as 



T (o) 



1 + AT (1) • I 



1 + + AS« • I 



1$,: 
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where T^ 1 ' and are the closed-shell and one-valence 
PRCC operators, respectively. The superscript (1) is 
used to indicate the perturbation. The diagrammatic 



representation of these cluster operators are shown in 
Fig. [2j 



value equation as 



Using Eq. (11) in Eq. (10), we can rewrite the eigen- 
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To derive the PRCC equations, we project the above 



equation on e 



and retain the terms linear in A. In ad- 



dition, for further simplification, we use normal-ordered 
form of the Hamiltonian if N = H - ($ V \H\$ V ). After 
these sequence of operations, the eigenvalue equation is 
modified to 



i? N S (1) + i? N T (1) (l + S (0) ) + Hl s c D (l + S (0) )] 1*0 
AE V S& + A£„T«(1 + |$03) 



where AE V = E v — (& V \H \<& v ), is the correlation en- 
ergy of the one- valence system. Like introduced ear- 
lier, H^ciP = e _T< ° H™ c e Tm is the similarity trans- 
formed NSD-PNC interaction Hamiltonian in the elec- 
tronic space. The PRCC equations of S 1 - 1 ' can now be 



derived by projecting Eq. ( 13 ) with the excited determi- 
nants and as 



NSD 
elec 
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While deriving the equations we have used the relations 
($P|TW|$„) = and (^T^S^O = 0. These follows 
as TW is an operator of closed-shell sector, it does not 
contribute to the PRCC equation of and S^. The 
closed-shell operators are the solutions of the similar 
set of coupled equations [35] 



(*g|{H N T( 1 )}|$ ) 
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- A^ T (1) |$ )(15b) 



These equations can be derived from the closed-shell per- 
turbed eigenvalue equation. We can also derive a simi- 
lar set of PRCC equations for the NSI-PNC interaction 
Hamiltonian. One major difference is, the cluster opera- 
tors are rank zero operators. 

After solving the RCC and PRCC equations, we can 
use the atomic states for the properties calculations. The 
RCC expressions and the diagrams contributing to the 
hyperfine structure (HFS) constants and the El transi- 



tion amplitudes are derived and discussed in our previous 
work [34] . In the present work, we use the same expres- 
sions and diagrams to compute HFS constants and El 
transition amplitudes. The PNC induced electric dipole 
transition amplitude, using PRCC wave function, is 



E1PNC = (*,„||D||*, 



(16) 



where D is the dipole operator. This expression, unlike 
the conventional sum-over-sates approach, implicitly ac- 
count for all the possible intermediate states. From Eq. 
( 11 1, for the NSD-PNC interaction, the transition ampli- 



tude is 
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Consider terms linear in A and retain only those up to 
second order in cluster amplitude. Define the electronic 
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component as £7l£{^?, corresponding to the H^JjP, ^ 1S 
then given as 
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($„,|DT« +T'°> t DTW + T« T DTW 

+T^T> + DT< W + T (1 ) f 5 (0)t D 
+5 (0)t DT (i) + T (i)t D5 (o) + DS (i) + g(i)t D 

+5 (o)t DS (i) +s (i)t D5 (o)| $i; ^ (18) 



To calculate E1PNC, we use diagrammatic analysis to 
identify the Goldstone diagrams from these terms. How- 
ever, we exclude the structural radiation diagrams, aris- 
ing from the terms involving two-body cluster operators, 
for example, T^DT 2 (0) . The selected diagrams from the 
leading order and next to leading order terms are shown 
in the Fig. [3j 
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FIG. 3. Some of the leading order PRCC diagrams which 
contribute to the _B1™ C C of one-valence atoms. 



III. RESULTS AND DISCUSSIONS 



A. Single-particle basis functions 



the spinor-spherical harmonics. The radial components 
are then defined as linear combination of Gaussian type 
functions [5rJl [57] 



p 

Q„«(r)=£C*4>). 



(21) 



The index p = 1, 2, . . . , to, where to is the number of basis 
functions and C' Kp are the coefficients of linear combina- 
tion. For large component we choose 



dL(r) = C l v 



(22) 



where n K is an integer and . is the normalization 
constant. The small component are derived from the 
large components using kinetic balance condition. The 
a p follow the general relation 



p-i 



(23) 



The parameters, ao and /3, are optimized such that the 
single particle energies of the core and valence orbitals are 
in good agreement with the numerical results, obtained 
from GRASP92 [35]. In Table. |TJ we compare the energy 
of the valence orbitals from the GTO with the GRASP92 
data. 



TABLE I. The valence orbital and SCF energies of Gaussian 
type orbitals (GTO) are compared with the GRASP92 data. 



Orbitals 



GTO 



GRASP92 



6s *S 1/2 
6p 2 Pi/2 
Qp 2 P 3 /2 
5d 2 D 3/2 
5d 2 D 5/2 
Escf 



-0.413668 
-0.301112 
-0.288305 
-0.303070 
-0.300885 
-14067.0622 



-0.413665 
-0.301113 
-0.288307 
-0.303071 
-0.300886 
-14067.0676 



For all the calculations we use Gaussian type orbitals 
(GTOs) or single particle wave functions with V N ~ 2 po- 
tential. As mentioned earlier, to incorporate the rela- 
tivistic effects we use the Dirac-Coulomb atomic Hamil- 
tonian. For the nuclear potential we consider the finite 
size Fermi density distribution 

ft^(r) = 1 ~ e ff_ c)/a , (19) 

here, a = t41n3. The parameter c is the half-charge 
radius, that is p n uc(c) = jOo/2 and t is the skin thickness. 
The orbitals are of the form 

ib (r) = -( P ««( r )x«"( r / r ) ^ (20) 

where P nK (r) and Q nK {r) are the large and small com- 
ponent radial wave functions, k is the relativistic total 
angular momentum quantum number and XKm{ r / r ) are 



B. Excitation energies, hyperfine structure 
constants and El transition amplitudes 

The excitation energies, hyperfine structure constants 
and the El transition amplitudes from our calculations 
are listed in the Tables. [TTJ |III| and |IV[ respectively. 
These results are obtained using a fairly large basis of 
177 active GTOs, it consists of 19, 17, 17, 17, 15 and 13 
orbitals in the s, p, d, f, g and h symmetries, respec- 
tively. To arrive at this basis set we start with a moder- 
ate size of 100 active orbitals with the combination 12s, 
lOp, lOd, 10/, 8g and 6h. And perform seven sets of 
calculations by adding one orbital to each symmetry in 
every successive sets. The % change in HFS constants 
and El transition amplitudes with respect to number of 
active orbitals are shown in Fig. [4] As we see in the 



figure, the El transition amplitudes converge and there 
is no observable change in the amplitudes after 155. On 
the other hand, for HFS constants we observe a slower 
convergence pattern. It is evident from the figure that 
the HFS results are close to convergence. The maximum 
uncertainty is about 0.5%, in the case of 5d 2 D 5 / 2l but 
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it is smaller for the states 6s 2 5 , 1 / 2 , 6p 2 P\/ 2 , 6p 2 P-. 
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FIG. 4. The convergence (% change) of the hyperfine struc- 
ture constants and the El transition amplitudes with respect 
to the number of active orbitals. 



The excitation energies from our calculations are listed 
in Table. [TTJ As described in Sec. [TTJ these are calculated 
using RCC with the CCSD approximation. Except for 
the 6p 2 P 3 / 2 excitation energy, our results are better or 
on par with the previous theoretical results when com- 
pared with the experimental data. The all-order results 
reported in Ref. [29] are closer to the experimental data 
than the other theoretical results, including the present 
work. For the 6p 2 P 3 / 2 , our result is in close to the 
RCC result of Sahoo and collaborators [57j. Among the 
other three results for this level, our result is closer to the 
Ref. [2§]. The dominant excitations which contribute to 
the denominator of the NSI-PNC matrix are 6s 2 5 , 1 / 2 — 
6p 2 Pi /2 and 6p 2 P 1/2 - 5d 2 D 3/2 . However for the NSD- 
PNC, these are 6s 2 S 1/2 - 6p 2 Pi/ 2 , Qp 2 Pi/2 - 5d 2 D 3/2 , 



6s 2 S 1/2 - 6p 2 P 3/2 and 6p 2 P, 



3/2 



3/2 



5d 2 D 3/2 . The ac- 



curacy achieved for these in the present work are 3.4%, 
4.5%, 2.4% and 4.6%, respectively. We have incorpo- 
rated these errors in the total uncertainty estimates for 
the PNC results. 

In the Table. |III| we present, and compare HFS con- 
stants obtained from the present calculations with the 
other theoretical and experimental results. As evident 
from the table, our results for the 6p 2 Pi/ 2 , 5d 2 D 3 / 2 
and 5d 2 D 5 / 2 states are in better agreement with the ex- 
perimental data than the other theoretical results. How- 
ever, for the 6s 2 5 , 1 / 2 state, like the other theoretical 



TABLE II. Excitation energy for some of the low lying exci- 
tations in 171 Yb + . The values are in cm" 1 . 



Level This work 



Other works 



Exp.RefgS]. 



3/2 



and 5d 2 D 3 / 2 , the uncertainties are 0.3%, 0.3%, 0.1% and 5^ 
0.05%, respectively. 



5/2 



Qp 



1/2 



6f>3/2 



23983 



25576 



27985 



31757 



23926 a 
21238 b 
22711 c 
22820 d 
22449 b 
24178 c 
24261 d 
28749 a 
28048 b 
27945 c 
27945 d 

28109(1000)° 
32376 a 
31411 b 
31403 c 
31481 d 

31604(800)° 



22961 



24333 



27062 



30392 



a Reference [28]. 
b Reference [40], 

c Reference [29 -MBPT + corrections. 
d Reference [29]-All-order. 
c Reference [27]. 



TABLE III. Magnetic dipole hyperfine structure constants of 
m Yb+ in the unit MHz. 



State This work 



Other works 



Exp 



6s 1/2 13488.314 



6p 1/2 2348.036 



6f>3/2 



5rf. 

5rf; 



3/2 
/2 



313.522 



421.131 



.567 



a Reference^. b 
d Reference [27] . e 
g Reference 02], h 



13217 a , 13172 b , 12645(2) e 
13091 c ,13332(1000) d , 
12730(2) e 

2533 a ,2350 b , 2104.9(1.3) c 
2371 c ,2516(400) d , 
2317 c 

388 a ,311.5 b ,330 c , 877(20) f 
322(20) d ,391 e 

291 a , 489 c , 447(20) d , 400. 5 g , 430(43) h 

-96 c ,-48(15) d , -12.6 s -63.6(7) s 

Reference 00], c Reference^. 
Reference [IT]. f Reference 03] . 
Reference [44]. 1 Reference [45 j. 



results, our result is also larger than the experimental 
result. Among the theoretical results, our result for this 
state is closer to the Ref. [27 . The reason for this is 
the method employed and the type of the orbitals used 
in the two works are similar. For the 6p 2 P 3 / 2 state there 
is a large discrepancy between the theoretical results and 
experimental data. However, it must be emphasized that 
the experimental data is from a relatively old measure- 
ment. Our result lies between the third-order MBPT 



G 



results of Ref. gD] and the RCC result of Ref. [27]. 

The impact of the electron correlation effects is dis- 
cernible in the Table. [V] where we list the contributions 
from various RCC terms. The RCC terms in the table 
are based on the expression in our previous work |34j . 
The contribution listed as "Other" correspond to the 
terms S , |i?hf s T + c.c. and S^H^TSi +c.c. As expected, 
the dominant contribution is from the DF term. It con- 
tributes approximately about 72%, 66%, 58%, 69% and 
110% for the states 6s 2 S 1/2 , 6p 2 P 1/2 , 6p 2 P 3/2 , 5c? 2 D 3/2 
and 5c? 2 D 5 / 2 , respectively. Our DF value 9716.7 and 
1548.2 for states 6s 2 Si/ 2 and 6p 2 Pi/ 2 , respectively are 
on the higher side of the values, 9577 and 1542, reported 
by Safronova and collaborators in their recent work [29"] . 
On the other hand, for the 6p 2 P3/ 2 , 5c? 2 D 3 / 2 and 
56? 2 L> 5/2 states our results of 182.5, 289.7 and 110.2, 
show a close match with the values 183, 290 and 111 
from Ref. [25] . The next two leading order contributions 
are from the terms Sj^hfe+cc. and S^-ffhfs + c.c.. Unlike 
other states, 5g? 2 D b / 2 shows a different correlation pat- 
tern and contribution from S 2 Hm s +c.c. is about -321% of 
the total value. However S'jiThfs+c-C. contributes 44% of 
the total value. Despite the large cancellations our total 
result compares well with the experiment. 



TABLE IV. The electric dipole transition amplitudes of 
171 Yb+. 



Transition 


This work 


Other works 


Exp. 


6pi/ 2 <- 


— 6S!/ 2 


2.748 


2.72 a ,2.73 b , 
2.75 c , 2. 64 d , 2.72(1)° 


2.47(3) 


6P3/2 <- 


— 6Sl/ 2 


3.901 


3.84 a ,3.84 b , 
3.83 c ,3.71 d ,3.83(l) c 


3.36(3) 


5g? 3 /2 <" 


~ Spi/2 


3.138 


3.09 a ,3.78 b , 


2.97(4) 








3.06 c ,2.98 d ,3.06(2) e 

1.36 a , 1.55 b 
1.35 c ,1.32 d ,1.35(2) e 




5g? 3 /2 «- 


— 6p 3/2 


1.369 




5c?5/2 *~ 


— 6p 3/2 


4.307 


4.77 b ,4.23 c 






4.23(3) e 





a Reference [28], 
b Reference [40], 

c Reference^- MBPT + corrections. 
d Reference 29 - All-order. 
e Reference [27]. 
f Reference [46] [47]. 
s Reference |48]. 



The El transition amplitudes are presented in the Ta- 
ble. IV For comparison the results from the other theo- 
retical and experimental works are also listed. Like HFS 



constant, transition amplitudes are calculated using the 
RCC wave functions. We have used the similar expres- 
sions and diagrams as the HFS except for one key dif- 
ference, the hyperfine operator is replaced by the dipole 
operator. The experimental results are available only 



6p 2 P\/ 2 — 5c? 2 I?3/2 transitions. Among all the theoreti- 
cal results, the results from the recent all-order work [29] 
are closest to the experimental data. All other results, 
including ours, are on the higher side of the experimen- 
tal value. The component wise contributions are listed in 
the Table. [V] Like in the case of HFS constant, the DF 
term has the dominant contribution. It contributes ap- 
proximately about 118%, 116%, 123%, 124% and 121%, 
respectively, for the transitions listed in the table. A 
close agreement is observed in the DF data from our cal- 
culation with the Ref. 1291. 



C. NSI-E1PNC 

For calculation of NSI-E1PNC, we use the expression 
Eq. (fl8l derived for NSD-E1PNC in the electronic space. 



However, the important difference in this case is, as men- 
tioned earlier, the PRCC cluster operators are rank zero 
operators. In terms of diagrams, the ones with dominant 
contributions are derived from Fig. [3] with the NSD- 
perturbed operators replaced by the NSI-perturbed ones. 
In Table. |VI[ we list the contributions from the different 
terms in PRCC. Among all the terms, the largest con- 
tribution, about 117% of the total value, is from DS^ . 
The reason for this, as evident from Table. |VII[ is the 
large £?pnc mixing between 6s 2 Si/ 2 and np "P\/ 2 or- 

bitals. This large contribution from DS^ is consistent 
with the pattern of correlation reported in Ref. [27] . The 
next leading order contributions are DS { 2 ] + H.c. and 
The former involve one core and one virtual or- 



(i) 



bitals, and the later connects T{ 
orbital. These contribute about 



tively. The terms DS[ 0) ■ 



and D through a core 
17% and 15%, respec- 



TI.c. and S{ 0)i DS{ 0) +c.c. 



are third and fourth leading order terms, contributing 
about -7% each. The contribution from normalization is - 
2.8%. Small but not insignificant contributions of 2% and 



-1.8% are also observed from the terms T 2 °^ ' DT[ L ' +c.c. 

and S^ D, respectively. 

To examine the correlation pattern more closely we 



and T^D, 



and for 



pick the leading order terms DS{ ^ 1 
these we calculate the E1PNC contributions from various 
intermediate np 1 / 2 states. The dominant contributions 
from these are tabulated in Table. IVIII The same anal- 

As 



VIII 



ysis but at the DF level is presented in Table, 
we see in both the tables, dominant contribution is from 
the 6p 2 P\/ 2 state, contributing about 117% of the total 
value. The reason for this is, large -ffpNC induced mixing 
with the energe tically closer 6s 2 S 1/2 state. The PRCC 
value is about 42% larger than the Dirac-Fock contri- 
bution. This can be attributed to the large amplitude 
of Si , and hence to the correlation effects incorporated 
using PRCC. The next dominant contribution among the 



for the 6s 2 S 1 / 2 — d>P 2 P\j2i 6s 2 S 1 / 2 — 6p 2 Py,/ 2 and core orbital is 5pi/ 2 . This contributes to Tj D through 
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TABLE V. Magnetic dipole hyperfine constants and El transition amplitudes, contributions from different terms in the RCC. 
The operator "O" here represents the hyperfine interaction Hamiltonian i/hfs for HFS constant and the dipole operator D for 
the El transition amplitudes. 

State/Transition Coupled-cluster terms 



DF O - DF S{6 S\d SlOSi SjOSi S ] 2 dS 2 Other Norm 

+C.C +C.C +CC. 



6s 1/2 9716.682 -427.318 2756.071 1145.559 125.344 204.287 243.111 -49.371 -225.956 

6p 1/2 1548.208 -49.951 528.950 245.862 28.191 46.824 23.105 22.266 -45.402 

6p 3/2 182.531 -5.430 56.715 53.520 5.802 4.570 18.932 1.897 -5.012 

5d 3/2 289.667 7.017 82.212 4.875 3.863 5.924 30.457 4.316 -7.197 

5d 5/2 110.234 4.187 29.781 -220.089 -16.443 2.032 19.468 1.186 1.076 



< — 6Si/ 2 


3.242 


0.001 


-0.175 


-0.311 


-0.010 


0.019 


0.029 


0.006 


-0.052 


< — 6Si/ 2 


4.543 


0.004 


-0.254 


-0.378 


-0.012 


0.023 


0.037 


0.006 


-0.067 


< — 6j>l/2 


3.861 


0.005 


-0.437 


-0.287 


-0.003 


0.032 


0.034 


-0.004 


-0.068 


< — 6p 3/2 


1.697 


0.002 


-0.206 


-0.121 


-0.000 


0.012 


0.013 


-0.001 


-0.027 


< 6p 3 / 2 


5.200 


0.010 


-0.566 


-0.326 


-0.001 


0.034 


0.036 


0.005 


-0.079 



the -ffpNC perturbed 6s 2 5 l 1/ / 2 - In this case as well PRCC 
contribution is larger than the DF. 

The total NSI-E1PNC result from our calculation is 
presented in Table. [X] We have also listed the DF con- 
tribution. The other two theoretical results are based on 
the calculations with correlation-potential-method [55J 
and RCCSD(T) [27]. The E1PNC result from these two 
works differ from each other substantially. The CCSD(T) 
result from Ref. [57] is about 26% larger than Ref. [28] . 
Our DF value is marginally on the higher side of the 
value reported in Ref. [27] . However, the total result lies 
between Refs. [55J and [57J, but closer to the coupled- 
cluster result of Ref. [57J . 



D. NSD-E1PNC 

For the NSD-PNC, the dominant contributions from 
various PRCC terms in Eq. (181 are listed in Table. 
|VI[ For hyperfine transitions F v = — » F w = 1 and 
F v = 1 ->• F w = 2, like in NSI-PNC, DSf 3 is the lead- 
ing order term. It contributes about 231% and -130%, 
respectively. For transition F v = 1 — > F w = 1, how- 



spectively. Non-negligible contributions are also observed 



ever, S^\D is the dominant term, contributing about 
-123% of the total value. The same trend is reported 
in Ref. [35], where the contributions are about 252%, - 
226% and -157%, respectively, for the F v = ->■ F w = 1, 
F v — 1 — > F w — 1 and F v = 1 —> F w = 2 transitions. 

The next leading order term, S^\d, contribute about 
-114% and 35%, respectively to the F v = F w = 1 
and F v = 1 — > F w = 2 transitions. However, the second 
leading order term is DS[^ for the F v = 1 — > F w = 1 
transition, it contributes about 69%. The next two lead- 
ing order terms are T^\d and DS^ +H.c The contri- 
butions from these terms, in the sequence listed in Table. 
VI] are 39%, 11% and -22%, and -24%, -59% and -4%, re- 



from the terms T 9 (0) ' DT^> + c.c 



.(0) 



S^'DS?' +c.c. 

Unlike NSI-PNC, in NSD-PNC np 3/2 states also con- 
tribute to the E1PNC matrix element. In Table. [IX] we 
list dominant contributions from odd parity r2.p1/2 and 
np 3 / 2 states in the PRCC calculations for the F v = — » 
F w = 1 transition, and specifically the contributions from 

DS^, SW[d and T^\d. At DF level we pr esent the 

As 



and 



(0) 



VIII 



contributions from the npi/2 states in Table, 
we see in these tables, both at DF and PRCU levels, 
the dominant contribution is from the 6p 2 P\/2 state. 
The total contribution from this in the PRCC calcula- 
tions is about 150%, which can be attributed to the 230% 



and -79% contributions from DS[ L) and S^\d, respec- 
tively. The large contribution through DS^ is due to 
the strong -ffpNC mixing with 6s 2 S1/2 ■ However that is 
not the case with 5d 2 D 3 / 2 , which contributes through 

S^\D, At DF level 6p 2 P\/2 contributes only through 
mixing with 6s 2 5*1/2, and contribution is about 214%. 
The state 6p 2 P 3 / 2 is the third most dominant contribut- 
ing one, contributing about -40%, through -ff PNC per- 
turbed 5d 2 D 3 / 2 - The other higher energy orbitals have 
negligible contribution. 

The NSD-E1PNC total results are given in Table. [X] 
For comparison we have listed the DF contributions. Our 
DF results 6.915, 1.632 and —3.643 for the three hyper- 
fine transitions listed in Table. [X] compare well with 
the results 6.90, 1.70 and -3.70, reported in Ref. [29]. 
Our total results, however, are on the higher side of the 
random-phase approximation (RPA) based results from 
Ref. 29 for all hyperfine transitions. The other theo- 
retical NSD-PNC data available for comparison is from 
Dzuba et al [55J using the correlation-potential-method. 
Our results for transitions F v = — > F w = 1 and 



TABLE VI. The NSI and NSD E1PNC component wise contribution from various terms in the PRCC. The NSI and NSD 
contributions are listed in the units of ieag x 1CP 11 (— Qw/N) and ieao^'w x 10 -12 , respectively. 



Transition 


DS^ 




DTP 


T m\ D 


DSP 
+cx. 


T (o)t OTi (i) 

+C.C. 


T (0)t OT (i) 

+CC. 


+C.C. 


+C.C. 


+C.C. 


Other 


Norm 


NSI-PNC 




8.950 


-0.139 


-0.005 


1.179 


— 1.278 


-0.025 


0.168 


-0.511 


-0.529 


0.028 


0.003 


-0.215 


F F 












NSD-PNC 














1 

1 1 

2 1 


6.689 
1.431 
-3.590 


-3.301 
-2.539 
0.952 


0.030 
-0.002 
-0.020 


1.118 
0.221 
-0.608 


— 0.683 
-1.221 
-0.114 


0.001 
0.000 
0.000 


-0.250 
-0.061 
0.130 


-0.239 
0.026 
0.163 


-0.279 
0.038 
0.194 


-0.008 
0.008 
0.009 


-0.099 
-0.021 
0.053 


-0.083 
0.059 
0.079 



TABLE VII. The NSI E1PNC dominant contribution from 
the intermediate odd parity states in the PRCC. The listed 
E1PNC values are in the units of ieao X 10 _11 (— Qw/N). 



Ea 



1/2 



E, 



6pi/2' 



The uncertainty associated with H] 



PNC 





DS 


(i) 
i 








Id 




D 


S$ x > 


E1PNC 


state 


D 


T (i)t 


E1PNC 


state 


-3.861 


100.787 


8.918 


GPl/2 


0.003 


-1.184 


0.0 


2P1/2 


-0.217 


-29.596 


-0.147 


7P1/2 


-0.010 


2.753 


-0.001 


3P1/2 


0.047 


16.932 


-0.018 


8P1/2 


-0.008 


7.623 


-0.002 


#1/2 


-0.009 


-21.986 


-0.005 


9P1/2 


1.290 


39.984 


1.182 


5P1/2 


0.106 


-29.329 


0.071 


10pi/2 








-0.161 


23.995 


0.088 


HPl/2 










-0.096 


15.017 


0.033 


12pi/2 











F v — 1 — > F w — 2 are in good agreement with their 
results. For transition F v — 1 — > F w — 1, however, our 
result is higher than their value. 



E. Uncertainty estimates 

To calculate the uncertainty of our E1PNC results 
we resort to a analysis based on the sum-over-state ap- 
proach. In this method, the net uncertainty associated 
with an intermediate state is 



A = 5E C 



8E1 + 5H } 



PNC, 



(24) 



where SE CKC1 and SE1 are the deviations of excitation en- 
ergy and El matrix element form the experimental data. 
And these are calculated based on results presented in 
Tables. |n|and IV respectively. For SHp^c, uncertainty 
associated with -ffpNC matrix, we resort to the deviation 
oi -^J AiAf from experimental data, where A t and Af rep- 
resent the magnetic dipole hyperfine constants of initial 
and final states of the E1PNC transition. 

As discussed earlier, dominant contribution to 

P\/2 state, which con- 

-f/pNC matrix el- 
^1/2)1 El matrix 
1 2) and energy denominator 



NSI-PNC is from the 6p 

tributes through DS[^ 
ement (6p 2 P 1 / 9 |i?i 



i.e. 

iwc|6s 



(5d 2 D 3/2 \D\6p 2 A 



matrix element, we get using our RCC results for hy- 
perfine constants, is 9.1%. And, the relative uncertainty 
of El matrix element and the energy denominator are 
calculated as 5.7% and 3.4%, respectively. Combining 
these, the net uncertainty in NSI-PNC result is 18.2%. 

For NSD-PNC as well 6p 2 Pi/2 is the dominant con- 
tributing state. However in this case unlike NSI-PNC, 

apart from DS^\ contribution through S^\d is not 
negligible. The matrix elements involve in this case are 
(5d 2 £ 3 / 2 |#PNc|6p 2 Pi/ 2 ) and (6p 2 Pi /2 \D\6s 2 5 1/2 ), 
and the energy denominator is Eq Pi/2 — E^d 3/2 . Using 
the similar analysis we get 4.52%, 11.26% and 2.41%, re- 
spectively for (Ji/piMCi SE1 and 5E. Combining these, we 
get 18.19% as the net uncertainty associated with term 

S^\D, The rms relative uncertainty in NSD-PNC re- 
sults are then 18.17%. 



IV. CONCLUSIONS 

In this work, we present NSI and NSD-PNC transition 
amplitudes for the [4/ 14 ]6s 2 S 1/2 - [4/ 14 ]5d 2 D 3/2 tran- 
sition in 171 Yb + ion. To estimate the uncertainty of the 
the PNC results, we also calculate excitation energies, hy- 
perfine structure constants and El transition amplitudes 
for some of the important low lying states, using RCC 
theory. The E1PNC results are computed using PRCC 
theory, which is formulated based on RCC theory and 
it incorporates electron correlation effects arising from 
a class of diagrams to all order in the presence of PNC 
interaction as a perturbation. 

Our results for excitation energies, hyperfine structure 
constants and El transition amplitudes are in good agree- 
ment, in some cases better, with the previous experimen- 
tal data. Our NSI-PNC result lies between the results 
from two previous studies reported in Refs. [55] and |27j . 
The NSD-PNC DF results from our work is in excellent 
agreement with the results reported in Ref. [29] for all 
hyperfine transitions. The total NSD-PNC result for the 
F v = — > F w — 1 hyperfine transition lies between the 
results of Refs. [25] and [2S]. For the remaining two, our 
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TABLE VIII. The Dirac-Fock dominant contributions from the intermediate odd parity states. The listed NSI-E1PNC and 
NSD-E1PNC values are in the units of ieao x 10 -11 (— Qw/N) and ieao x 10~ 12 (i' w , respectively. 













i/pNC-D 






D 


-ffpNC 


E1PNC 


state 


D 


-ffpNC 


E1PNC 


state 


NSI-PNC 


-3.861 
-0.217 
0.047 
-0.009 


71.073 
-18.657 

10.547 
-13.620 


6.288 
-0.092 
-0.011 
-0.003 


6pi/ 2 

7pi/2 

8pi/ 2 

9P1/2 

NSD-PNC 


0.003 
-0.010 
-0.008 

1.290 


-1.178 
2.657 
6.738 

23.438 


0.0 
-0.001 
-0.001 
0.693 


2pi/2 
3pi/ 2 

4pi/2 

5pi/ 2 


-3.861 
-0.217 
0.047 
-0.009 


-118.154 
31.016 

-17.534 
22.642 


6.210 
-0.092 
-0.011 
-0.003 


6pi/2 
7P1/2 
8pi/2 
9pi/2 


0.003 
-0.010 
-0.008 

1.290 


1.957 
-4.417 
-11.202 
-38.965 


0.0 
-0.001 
-0.001 
0.684 


2pi/2 

3?>l/2 
4pi/2 
5P1/2 



TABLE IX. The NSD E1PNC dominant contribution from the intermediate odd parity states in the PRCC. The E1PNC values 
are in the units of ieao x I0~ 12 fi' w . The contributions listed from the core riPx/i an d nP 3 / 2 orbitals are from the terms T m [D 
and DT^\ respectively. 





DSP 










Orbital 




T m \D/ DT W 




Orbital 


D 




E1PNC 


D 




E1PNC 




D 


T (i)t 


E1PNC 




-3.861 
-0.217 

0.047 
-0.009 

0.106 
-0.161 
-0.096 


-126.570 
37.839 

-21.776 
28.408 
38.534 

-32.883 

-22.179 


6.653 
-0.112 
-0.013 
-0.004 
0.056 
0.072 
0.029 


3.242 
-093 
0.011 
0.013 
0.075 
-0.081 
-0.053 


78.106 
-9.111 
4.704 
-5.743 
-6.647 
4.431 
1.876 


-2.298 
-0.008 
-0.001 
0.001 
0.005 
0.003 
0.001 


6P1/2 

7pi/2 
8pi/2 
9P1/2 
10pi/2 
llPl/2 

12pi/ 2 


0.003 
-0.010 
-0.008 

1.290 


1.965 
-4.550 
-12.416 
-63.798 


0.000 
-0.001 
-0.001 

1.120 


2P1/2 

3pi/ 2 
4pi/ 2 

5P1/2 



1.697 


-8.244 


0.000 


-4.543 


30.184 


-0.984 


0.024 


2.418 


0.000 


0.358 


-5.101 


-0.013 


0.008 


-1.376 


0.000 


-0.140 


2.786 


-0.003 


-0.028 


1.845 


0.000 


0.136 


-3.670 


-0.004 


0.075 


-2.226 


0.000 


-0.050 


4.445 


-0.002 


0.080 


-1.527 


0.000 


0.022 


3.370 


0.001 


0.043 


-0.791 


0.000 


0.035 


1.675 


0.000 



6P3/2 
7P3/2 
8P3/2 
9P3/2 
10p 3 /2 
llP3/2 

12p 3 /2 



-0.001 

0.006 
0.046 
0.749 



-0.001 
0.002 
0.109 

-3.964 



0.000 
0.000 
0.000 
0.021 



2p 3 /2 

3P3/2 
4p 3 /2 
5P3/2 



TABLE X. The total NSI (in the unit ieao x 10 -11 (-Qw/N)) and NSD (in the unit iea fl' w x 10" 12 ) E1PNC results compared 
with the previous theoretical results. 



Transition 






This work 




Other works 






DF 




PRCC 




NSI-PNC 


(5d 3 / 2 | <- (6s 1/2 | 




7.002 




7.626 


6.262(20) a ,8.470 b 


NSD-PNC 












{5d 3/2l F w = 1| <r- (6s 1/2 ,F v 


= 0| 


6.915 




2.896 


3.1(1.9) a ,2.6 c 


{5d 3/2 ,F w = lj (6^1/2,^ 


= 1| 


1.632 




-2.061 


-1.3(4) a ,-1.5 c 


{5d 3/2 ,F w = 2| •(- {6s 1/2 ,F v 


= 1 


-3.643 




-2.753 


-2.6(1.3) a ,-2.2 c 



a Reference [28], b Reference [27]. c Reference [29], 
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results are slightly on the higher side. The upper bound 
to the theoretical uncertainty associated with E1PNC re- 
sults is about 20%. 
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